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The zero point kinetic energy of protons in water is large on the scale of chemical interaction energies( 
29 Kj/mol in bulk room temperature water). Its value depends upon the structure of the hydrogen 
bond network, and can change as the network is confined or as water interacts with surfaces. These 
changes have been observed to be large on a chemical scale for water confined in carbon nanotubes 
and in the pores of xerogel, and may play a fundamental, and neglected, role in biological processes 
involving confined water. We measure the average momentum distribution of the protons in salmon 
Na-DNA using Deep Inelastic Neutron Scattering, for a weakly hydrated ( 6w/bp) and a dehydrated 
fiber sample. This permits the determination of the change in total kinetic energy of the system per 
water molecule removed from the DNA and placed in the bulk liquid. This energy is equal, within 
errors, to the measured enthalpy for the same process, demonstrating that changes in the zero point 
motion of the protons, arising from changes in structure as water molecules are incorporated in the 
DNA, are a significant factor in the energetics of the transition from the A to B phase with hydration, 
in this case, providing the entire binding energy of the water molecules to the DNA. The shape of 
the momentum distribution in the dehydrated phase is consistent with coherent derealization of 
some of the protons in a double well potential, with a separation of the wells of .2A. 

PACS numbers: 



Chemical interactions occurring in bulk water typically represent small changes in the energy of constituents com- 
pared to the energy sequestered in the zero point motion of the protons in the water, primarily in that of the stretch 
mode. To the extent that this energy does not change, it may be ignored, and usually is. Nearly all simulations of 
water in biological systems are done with models of water for which these changes cannot occur. The energy does 
change, however, as the structure of the hydrogen bond network changes. With the advent of Deep Inelastic Neutron 
Scattering measurements [? ], it has become possible to observe the changes in zero point kinetic energy directly. 
These changes can be largc(>KT)and in either direction. Water confined in carbon nanotubes has a transition in 
which the kinetic energy of the protons in each water molecule is reduced by 6Kj/mol[? ] , while water confined in the 
smallest pores of xerogel the kinetic energy increases by 14Kj/mol[? ]. We present here a system where these changes 
play a significant, indeed, the dominant role, in an important and well studied biological process, the transformation 
of DNA from the A to the B phase with the absorption of water. 

DNA is hydrophylic, changing its structure as it incorporates water molecules, transforming continuously from the 
A phase to the biologically active B phase as the number of hydrated water molecules increases from 1 to about 20. [? 
] The energy input necessary to remove a molecule of water at any level of hydration and place it in bulk water can 
be measured with Differential Scanning Calorimetry. It varies with the level of hydration and the composition of the 
DNA, and is on the order of lOkJ/mol/water molecule, decreasing with increasing numbers of hydrating waters. [? ] 
The attraction that leads to the hydrophylicity, and the source of the energy that produces the deformation of the 
A structure to become the B structure, is usually thought of as primarily electrostatic. [? ] However, as , as the 
configuration of the surroundings of the water molecule change, the potential energy surface of the proton changes. As 
a consequence, the ground state wave function changes, and with that there is a concomitant change in the zero point 
kinetic energy of the proton. This can be observed indirectly in the red shift of the vibrational modes from the gas 
phase to the bulk liquid, where this change in zero point energy accounts for 25% of the enthalpy difference between 
the two phases at room temperature. It has recently become possible to measure accurately the kinetic energy of 
protons in water as well as other hydrogen bonded condensed matter systems. [? ] This energy is almost entirely due 
to the zero point motion of the protons. The temperature plays a relatively small role, as the kinetic energy depends 
primarily on the degree of localization of the protons in the structure surrounding them. As the water of hydration is 
incorporated into DNA, the structure surrounding both the protons in the water and those in the DNA changes, and 
we would expect changes in the kinetic energy of the protons involved. We show here that reductions in the kinetic 
energy of the zero point motion of the protons in weakly hydrated salmon DNA, containing 6 water molecules/base 
pair, are not only significant, they account for the entirety of the enthalpy of hydration of the molecules, at that level 
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of hydration. Changes in the zero point energy of the protons in the DNA-water complex are therefore a major factor 
in the transformation of the A to the B phase of DNA, at least at these relatively low levels of hydration. There is 
also evidence, in the shape of the momentum distribution, that some of the hydrogen bonds in the dry DNA, have a 
double well character, as revealed by a characteristic oscillation in the momentum distribution. This property is lost 
upon hydration. 

The kinetic energy of the protons in the DNA-water complex is readily measured using Deep Inelastic Neutron 
Scattering(DINS). This method has been described in detail elsewhere. [? ] We briefly summarize it here. At high 
energy transfers, the usual neutron scattering function S(q, w), which gives the differential crossection for scattering, 
is described by the impulse approximation limit, in which the target particle behaves as a free particle for the duration 
of the scattering process. It is given by 



/hq p 
n(p)S(uj- — -q-—)dp (1) 



where n(p) is the momentum distribution of the protons, fko and Hq the transferred energy and momentum, respec- 
tively, and M the mass of the proton. S(q,w) is thus the Radon transform of n(p), and as the transform is invertible, 
n(p) can be obtained directly from the scattering data. This inversion is accomplished by a method [? ] in which the 
momentum distribution is represented, for samples such as ours that are isotropic , as a series expansion 



where L^(^-j) are the associated Laguerre polynomials, and the parameter a is related to the kinetic energy of the 

proton by K.E.=^-ct 2 .[? ] The coefficients a n are determined by a least squares fit of the data for S(q,u>) The 
instrument used was Vesuvio, at ISIS, the pulsed neutron source at the Rutherford Laboratory in England. The 
instrument has 64 detectors, arranged with scattering angles from 35 to 75 deg on both sides of the beam. Since in 
the impulse approximation limit, the data is determined entirely by the scaled variable y = ^ (u> — ^jj^j , all detectors 

give the same information when the time of flight data is converted to the y variable. The parameters of the fit are 
obtained by simultaneously fitting the data from all the detectors. The consistency of the data can be checked by 
comparison of the rescaled data from individual detectors, which all collapse onto a single curve. 

The measurements were made on fibrous salmon Na-DNA obtained from Sigma. The fibers were mechanically 
broken up into segments 2mm in size to insure a random orientation in the sample. The total sample, approximately 
5 grams, was weighed and packed into an aluminum container. No special attempt had been made to control the 
humidity of the sample environment, which had been the local atmosphere for some weeks. The container was sealed 
with an indium seal, and placed in the neutron beam in an evacuated chamber. After obtaining the desired number of 
counts, the sample was removed from the container, weighed in air, and placed in a vacuum at 40 degrees Centigrade 
for 24 hours. The sample was reweighed, again in air, and packed into the same container and resealed. The weighing 
process took approximately 20min. After counting for a day, the sample was removed and promptly reweighed. The 
difference between the weight upon removal from the dehydrating oven and upon removal from the sample holder was 
.012gr, presumably due to hydration during the period the sample was being packed into the sample holder or being 
weighed the second time, indicating that the protocol was sufficiently accurate as to the weight of the sample for our 
purposes. The weight of the initial hydrated sample was 4.76gr, that of the dehydrated sample as removed from the 
sample holder 4.16gr. We will use these values to calculate the number of water molecules/nucleotide removed by the 
dehydration process. This is accomplished using the known stoichiometry of the Na-Dna salt. We used the value for 
the ratio of G-C to A-T of .436. [? ? ] Our results are very insensitive to this ratio , as the weight ratio of protons to 
heavier atoms in the G-C and A-T pairs is nearly the same. The number of water molecules removed is 6/nucleotide, 
which is sufficient to produce a substantial change in the structure of the DNA. [? ] This is the number expected for 
the typical levels of relative humidity experienced by the sample. [? ] 

A comparison of the radial momentum distributions, 47rp 2 n(p), of the protons in DNA with (Hydrated) and without 
the hydrated water molecules(Dry) is shown in Fig 1. The values of a and the significant parameters in the fits are 
shown in Table 1 for the hydrated and dry samples. The error bars shown in the graphs are calculated from the 
least squares fitting procedure, and include the correlations between the parameters. Terms up to n=7 were included 
in the initial fits. The parameters set to zero were found to have average values less than half their rms variance, 
and did not significantly raise the chi-square of the fit when set to zero. Their inclusion enhances the rms point by 
point error, but has a negligible affect on the values of n(p).[? ] We can see immediately that the dry sample has 
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FIG. 1: Comparison of the average momentum distribution of the protons in hydrated Salmon DNA with that in dry Salmon 
DNA. The oscillation in the distribution of the dry sample indicates coherent derealization of some of the protons in a double 
well. The position of the local minimum allows an estimation of the separation of the wells as .2A 



TABLE I: The fitting coefficients and their variances for the momentum distribution measured in dry and hydrated DNA 



Parameters of measured n(p) 





Hydrated 


Dry 


a (A~ L ) 


4.81± .028 


5.19 ± .067 


n 


a n 


8a n 


a„ 


Sa n 


2 







.657 


.134 


3 


-.259 


.075 







4 







-.438 


.177 


5 












6 


.086 


.101 


.354 


.178 



a larger second moment, as is confirmed by the values of a in Table 1. Moreover, n(p) is not monotonic for large p. 
The appearance of an oscillation in the momentum distribution is characteristic of a population of protons that are 
coherently distributed in the ground state of a double well potential. The bimodal spatial wavefunction, when fourier 
transformed to give a momentum wave function, produces oscillations with a wavelegth in momentum space of 7r/d 
[? ], where d is the separation of the wells. In this case, we obtain a value of d of .2A 

Evidently, the hydration of DNA, at least for the first 6 molecules of water, is accompanied by a reduction of the 
zero point kinetic energy of the water-DNA system. That energy is available to produce a distortion of the A phase 
of the DNA towards the B phase. We wish to see how significant it is compared to the enthalpy change per water 
molecule upon removing the same number of water molecules from the DNA and placing them in bulk water. This 
has been measured by K. Virnik et al.[? ] and others [? ? ] using Differential Scanning Calorimetry.The value is 
3.0±.15 KJ/mol for salmon DNA at this level of hydration. To do this, we compare the total kinetic energy of the 
dry DNA together with 6 water molecules in bulk water, with that of the hydrated DNA. 

That is, if the number of protons in the dry sample is Ndna , and the number of protons in the water 
that was removed from the hydrated sample is N# 2 o, then in the hydrated sample the total kinetic energy is 
{^H 2 o J r^DNA)'^^'wef The kinetic energy of the dry sample is ^dna ffy &Dry to which must be added the ki- 
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netic energy of the water molecules in bulk water, ^h 2 oj^[< j h 2 o- The change in the kinetic energy/proton, averaged 
over all the protons, from the dry state to the hydrated state is therefore. 

3h 2 

= — [N DNA a 2 Dry + N H2 ov 2 H2 o - (Nona + N H2 o)c7 2 Wet )] (3) 
The change in the kinetic energy per proton is therefore 

AK.E/N Tot . = ^[^(rtry - + ^(°ko - - 2 wJ (4) 

where Nxot — Nona + Nh 2 o- The sigma parameter for bulk water was measured at the same time as the other 
experiments descrived here. The value of er=4.76±.035A is consistent with earlier measurements^? ]. With the ratio 
of protons in the dry DNA to that in the hydrated DNA as calculated from the stoichiometry as .696, and the values 
for the a parameters and their uncertainty given in Table 1, we find the average change in the kinetic energy/water 
molecule, to be, since there are two protons in a water molecule 

2AK.E/N Tot . = 3.02 ± MKJ/mol (5) 

That is, the change in kinetic energy of the protons in the combined water-Dna system, in going from the dry phase, 
with the water molecules in bulk water and the DNA unhydrated, to the hydrated phase, completely accounts for the 
measured enthalpy change. If there is some water still in what we have been calling the dry phase, the change in kinetic 
energy/proton would be increased abve the value given above. This does not mean that the electrostatic interactions 
are insignificant compared to the kinetic energy changes, only that the reduction in the electrostatic potential energy 
from the incorporation of the molecules of water in the DNA is nearly equal to the increase in the elastic potential 
energy of the A phase as it deforms into the B phase. The reduction in zero point energy then provides the binding 
energy for the water molecules. 

The extent to which this change in zero point energy can be attributed to the protons in the hydrating water, rather 
than the hydrogen bonds in the DNA is unknown. Comparable changes in the kinetic energy of protons in water 
molecules have been seen in water confined in carbon nanotubes[? ] and in xerogel[? ], so one might think most of 
the change is due to changes in the environment of the protons in the water molecules. The shape of the momentum 
distribution for the dry phase, however, suggests changes in the base pair hydrogen bonds, as the potentials for the 
protons there are sensitive to the separation of the base pairs, and this is easily changed with a structural deformation. 
It is also the case that the measured changes in the stretch and bending modes vibrational frequencies for the hydrated 
water, which account for most of the zero point energy, are far too small( at most 1.5%) to account for the changes in 
the kinetic energy if we assume the protons in the water remain in approximately harmonic potential wells. [? ? ] In 
any case, our results show that a correct description of the energetics of hydration of DNA requires that the quantum 
derealization of the protons in the system be accounted for. We expect this will be generically true of other systems 
in which the hydrogen bond network of water is strongly distorted, due for instance, to the confinement of the water 
in cells. 
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